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Abstract 

We compute the pion Generalized Parton Distribution (GPD) in a 
valence dressed quarks approach. We model the Mellin moments of 
the GPD using Ansatze for Green functions inspired by the numerical 
solutions of the Dyson-Schwinger Equations (DSE) and the Bethe- 
Salpeter Equation (BSE). Then, the GPD is reconstructed from its 
Mellin moment using the Double Distribution (DD) formalism. The 
agreement with available experimental data is very good. 

Keywords: GPD; pion; Dyson-Schwinger equations; Impulse approxi¬ 
mation; Soft Pion Theorem 


1 Introduction 

Introduced in the 1990s BEIE], GPDs have been intensively studied both 
theoretically and experimentally. Concerning the proton GPD, several mod¬ 
els have emerged a El El [71 El [9] based on different kind of parametrizations, 
and fitted to data. Here we focus on the pion GPD, which we model in 
an original way through the triangle approximation, but using propagators 
and vertices coming from the numerical solutions of the DSE and BSE. This 
approach have been successful in the case of the pion Parton Distribution 
Amplitude (PDA) [10]. In the first section, the details of the model will be 
given. Mellin moment will be computed using functional forms coming from 
the solutions of the DSE and BSE. The results will also be compared with 
available experimental data. In the second sections, the main point will be 
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the reconstruction of the GPD itself thanks to its Mellin moments. We will 
details the DD method which reveals itself very convenient. In the third 
section, symmetry questions will be highlighted allowing us to go beyond the 
triangle diagram approximation. 

2 Computing the Mellin moments 

In the pion case, the GPD H is formally defined as the Fourier transform of 
a non local matrix element: 


0,2x=0 

( 1 ) 

Instead of modeling directly the GPD H{x, t), we focus on the compu¬ 
tation of its Mellin moments defined as: 

^n(^, = y x^H{x, t) (2) 

and which are even polynomials in ^ of degree at most n -|- 1. In order to 
model those object, the so-called triangle diagram approximation is used 
here, as depicted on Fig. [H Gomputing this diagram naturally leads to an 
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Figure 1: Model of GPD Mellin moment in the triangle approximation for a 
valence dressed quaxk.Left panel', quark case. Right Panel, anti-quark case. 
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expression for the Mellin moments as: 


/ Mh A A A 

(k ■ nriT,(k - P - -) S(fc - -) 

S7-n S(k+ j) ir,(k + j,P+^) S(k- P),(3) 

where 


S{p) 

= [ - *7-p + M] Am(p^), 

(4) 

Am(s) 

1 

s + M2 ’ 

(5) 

P^(/c,p) 

= z75—M^^ / dzp^(z) ] 

JTV J -1 

(6) 

Pu{z) 

= R.{i-z^r, 

(7) 


have been inspired by the numerical solutions of the DSE-BSE[T0]. P) 
is the amputated Bethe-Salpeter pion vertex, S{k) is the quark propagator 
and Tr stands for the trace on color, flavour and Dirac indices. M is a free 
parameter which can be seen as the effective quark mass in this algebraic 
model, u is also a free parameter, but in all the following work, it will be 
hxed to 1, a value which allows one to recover the pion asymptotic DA [10]. /tt 
is the pion decay constant and Ry is a normalization constant. Comparison 
with experimental data can be done on the form factor mn which is the 
Mellin moment of order 0 of the GPD, and the PDF. Unfortunately, no data 
at 7^ 0 is available now, despite some pioneering work [13]. The agreement 
with experimental data shown on Fig. |2]is very encouraging for this algebraic 
model [H] . 



Figure 2: Left Panel. Comparison of our model with experimental data for 
the pion form factor. Right panel: Comparison of our model with an PDF 
£t to data, before and after evolution at low scale. 
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3 Double Distributions 

DDs have been introduced very early in order to model GPDs[I], |T5l |16] 
because they ensure the polynomiality property of the Mellin moments. We 
remind the reader that the link between the GPD H and the DD is a Radon 
transform: 

H{x,^,t)= f d(3da 5{x — (3 — a^) {F((3,a,t) + (8) 

J\a\+\0\<1 

It is possible within our model, to rewrite the Mellin moments of the GPD 
in terms of Mellin moments of the DDs. From there, one can identify for any 
order of the Mellin moments the DDs, and thus get back the GPD through 
the Radon transform [T3|. Gonsequently, one can get an analytic expression 
for the GPD, and thus check the fundamental properties as polynomiality 
of the Mellin moment, continuity at x = ^ and the support property. All of 
them are fulfilled in our approach as it can be noticed on Fig. |3] 



l-x 


Figure 3: Left Panel Full GPD reconstruction at vanishing t.This plot illus¬ 
trates the properties of support and continuity. Right Panel. Behaviour of 
the PDF at large x. We get back here the pertubative behaviour. 


4 Symmetry considerations 

4.1 PDF case 

Getting an analytic expression for the GPD allows also one to check the 
consistency of the different assumptions. One of them is the fact that due 
to a combination of isospin symmetry, charge symmetry and the two-body 
nature of our problem (two dressed quarks), the PDF must be symmetric 
with respect to x = |. As it can be seen on left panel of Fig. 01 this 
symmety property is slightly broken. This due to the assumption to work 
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Figure 4: Left Panel. PDF obtained with the triangle diagram only. The 
symmetry with respect to x = | is slightly broken. Right Panel. Additional 
contribution beyond the triangle diagram. The sum of the two contributions 
is symmetric 


with a triangle diagram[T7]. It is also possible to take into account the 
twist two operators directly when computing the Bethe-Salpeter amputated 
vertex. Adding such a contribution gives back a symmetric PDF, as shown 
on Fig. m 

4.2 Soft pion theorem 

There is an interesting link between the pion GPD and the pion PDA [T8] . 
Denoting the pion PDA by ^(z), one can show that: 

= = = (9) 

Our model failed to reproduce this property. There is no mystery behind that. 
Indeed, the triangle diagram is consistent enough with the rainbow ladder 
truncation of the DSE, such that it will give automatically the soft pion theo¬ 
rem, provided that the Axial-Vector Ward Takahashi Identities (AVWTI) are 
fulhlled [19]. In this algebraic model, it is not the case. But it is enough to 
upgrade the model by using directly the numerical solutions of the BSE-DSE 
to fulhll the AVWTI, and thus the soft pion theorem. 


5 Conclusions and outlook 

Our study has shown that it is indeed possible to build a new model of 
GPD based on the Dyson-Schwinger equations which will respects most of 
the theoretical constraints. It has been shown that, it was possible to pursue 
all the computations analytically providing that one uses the DD formalism. 
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Furthermore the very satisfactory comparison with experimental data is a 
key encouraging point to go beyond this algebraic model and build another 
model using the numerical solutions of the DSE and BSE. It also pave the 
way for a proton GPD, based on a quark-diquark model. 
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